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Abstract 

We study the semilinear partial differential equation (PDE) associated with the non-linear BSDE 
characterizing buyer’s and seller’s XVA in a framework that allows for asymmetries in funding, 
repo and collateral rates, as well as for early contract termination due to counterparty credit risk. 
We show the existence of a unique classical solution to the PDE by first proving the existence 
and uniqueness of a viscosity solution and then its regularity. We use the uniqueness result to 
conduct a thorough numerical study illustrating how funding costs, repo rates, and counterparty 
risk contribute to determine the total valuation adjustment. 


1 Introduction 

We study the total valuation adjustment (XVA) of a European style claim under differential borrowing 
and lending rates. We adopt the framework introduced in the companion paper of lichuch et al. 
(2011 ), where the authors define buyer’s and seller’s XVA in terms of solutions to non-linear BSDEs. 
These characterize the portfolio process replicating long and short positions in the claim, and take 
into account early termination of the contract due to counterparty credit risk, funding spreads, as 
well as collateral servicing costs. When funding spreads are zero, security lending and borrowing rates 
coincide, and collateral rates paid by the taker equal the ones received by the provider, Bichuch et 
(2015) give explicit expressions for the XVA as well as for the corresponding replicating strategies. 
However, in most realistic market scenarios, rates are asymmetric and funding spreads constitute a 
significant driver of XVA. Despite the unavailability of closed form expressions, in this paper we show 
that we can exploit the connection between the BSDEs and the corresponding nonlinear PDEs to 
study how funding costs, collateral and counterparty risk affect the total valuation adjustment. 

Our study extends previous literature along two directions. First, we develop a rigorous study of 
the semilinear PDEs associated with the nonlinear BSDEs introduced in Bichuch et al. (201! ) and 
yielding the XVA prices. Previous studies in this direction include Jercurio (201 ), who gives the 
PDE representations for the lower and upper bound prices of options in an extended Piterbarg (2010) ’s 
model where borrowing and lending rates differ. However, he does not take into account counterparty 
risk. This is accounted for in Burgard and Kjaer (2011a) and Burgard and Kjaer (2011b), who consider 
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an extended Black-Scholes framework where two corporate bonds are introduced in order to hedge the 
default risk of the trader and of his counterparty. They generalize their framework in Surgard and 
Kjaer (2013) to include collateral mitigation and evaluate the impact of different funding strategies. 
We also mention Crepey and Song (201 ) who study a BSDE with random terminal time of the type 
appearing in a funding valuation context, but do not develop connections to the resulting PDE. 

We prove the existence and uniqueness of a viscosity solution to the PDE. We extend standard 
results of )elong (2013) by using the locality of the viscosity solution, and also by accounting for 
the possibility of default risk. We then show that the PDE also admits a classical solution implying 
existence and uniqueness of the classical solution. Our result contrasts with earlier studies in the 
funding literature, e.g. Burgard and Kjaer (2011 ), where smoothness of the price of the derivative 
contract with respect to the underlying spot asset price is usually assumed to hold, and the PDE is 
derived after assuming that the market does not admit arbitrage. Using the classical solution, we 
compute the replicating strategies of the traded claim in a portfolio consisting of a stock, and risky 
bonds of the trader and of his counterparty. To the best of our knowledge, this is the first study 
to establish this result. Crepey (2011 ) develops a connection between the funding BSDE and the 
solution of a classical PDE, but in a different framework where the resulting PDE is linear and the 
market model complete, see Section 5.2 therein. We also mention the recent work of lie and Rutkowski 
(2013) who establish classical solutions to quasi-linear PDEs corresponding with BSDEs yielding the 
fair bilateral prices of a European contingent claim. Their work, however, does not consider the 
possibility that the trader or his counterparty can default. 

A second noticeable contribution of our study is a comprehensive numerical analysis made possible 
by the previously established existence and uniqueness result. Earlier studies on funding, such as 
Bielecki and Rutkowski (2014), Crepey (2015a), and Brigo et al. (2012) lack the computational com¬ 
ponent. We find strong sensitivity of XVA to funding costs, counterparty risk and collateralization 
levels. Viewing both buyer’s and seller’s XVA as functions of collateral levels defines a no-arbitrage 
band whose width increases with the funding spread. As the position becomes more collateralized, 
the value of the closeout payment increases and leads the trader to implement a riskier strategy. This 
yields increased buyer’s and seller’s XVA. Interestingly, both of these quantities may decrease in the 
counterparty’s default intensity as the default premium compensation demanded by the investor for 
holding counterparty bonds becomes higher than the funding costs incurred for replicating the close¬ 
out position. Moreover, increased values of counterparty’s default intensity reduce both seller’s and 
buyer’s XVA by a similar amount. 

Our findings serve as a useful guide to risk-management and trading desks, who can decide on the 
terms of the trades (collateralization levels, borrowing rates charged by treasury, etc...) based on the 
incurred costs as measured by XVA. 

The rest of the paper is organized as follows. Section 2 reviews the framework for computing XVA 
introduced in iichuch et al (201 ). Section 3 analyzes the semilinear PDEs. Section 4 develops a 
comprehensive numerical analysis. Section 5 concludes the paper. 

2 The Setup 

We review the setup introduced in the companion paper by Bichuch ct al (2015) and refer to sections 
2 and 3 of their paper for additional details. We give here a self-contained exposition needed to 
understand the PDE analysis presented in the next section. We fix a probability space (P,t?,P), 
where P denotes the physical measure. The market filtration is Qt '■= V Wt, t > 0, where the 
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filtrations J~t and Tit will be specified later in the section. 

We take the viewpoint of an investor “I” (also referred to as trader or hedger) who wants to compute 
the total valuation adjustment, abbreviated as XVA, of a European claim traded with a counterparty 
“C”. The claim is written on a stock whose time t price is denoted by St- The claim matures at T and 
has terminal payoff where $ : R>o —> H is a piecewise continuously differentiable real valued 

function of polynomial growth. The value of the claim as determined by a third party evaluator is 
given by 


V(t,S t ) = e -^(T-t) E Q[$ (5T ) | jr t ] 5 

where the valuation measure Q associated with the publicly available discount rate rp chosen by the 
valuation agent is equivalent to P and given by the Radon-Nikodym density 
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2.1 The Portfolio securities 

This section describes the securities at disposal of the investor to construct his replicating portfolio. 
They include the default-free stock on which the financial claim is written, and two risky bonds 
underwritten by the trader and by his counterparty. Under the physical measure, the dynamics of the 
stock price process is given by 

dS t = nS t dt + a S t dWf , 

where /r and a are constants. Moreover, W v is a standard Brownian motion under the probability 
measure P. We set P := (Pt)t>o to be the (Q. P)-augmentation of the filtration generated by W v . 

Let rj and tc be the default times of the trader and of his counterparty, respectively. These are 
independent exponentially distributed random variables with constant intensities hj , i E {I, C}. We 
use Hf = l{ Ti< t}, t > 0, to denote the default indicator process of the i-th name. The default event 
filtration is set to Ht := <r(H\u <t). 

The two risky bonds underwritten by the trader I and by his counterparty C mature have the 
same maturity T of the claim. For 0 < t < T, i E {I, C}, their price processes admit dynamics given 
by 

dPl = ( r * + h?)Pi dt - P\_ dHl P* = e ~^ +h ^ T , (2) 

where r l is constant. Throughout the paper, we set r := 77 A tq A T. We also recall the following 
relation, given in lichuch et al (2015), between the physical measure P and the valuation measure Q: 

W t Q = Wf + — t 

(7 

™t Q = + J o (! - HI)(hi ~ hf)du, 

hf = h p + r l - r D , 

2.2 Security, Funding, and Collateral Accounts 

This section describes the various trading accounts used by the hedger to finance the strategy repli¬ 
cating the price process of the claim. 


ie{I,C}. (3) 
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2.2.1 The Security Account 

Financing and lending activities related to the stock security happen via the repo market. We denote 
by rj the rate received by the hedger when he lends money to the repo market. Symmetrically, we 
denote by r~ the rate he has to pay when he borrows money from the repo market. Let ipt be the 
number of shares of the security account at time t. The value of this account at t is given by 

Bj := Bj('ip) = efo ( 4 ) 

where 

r r {x) = r~t {x<0} + rjt {x>0} . (5) 

2.2.2 The Funding Account 

The hedger finances the portfolio strategy to replicate the claim from his treasury. We denote by rj 
the rate at which the hedger lends cash to the treasury, and by rj the rate at which he borrows cash 
from it. Let denote the number of shares of the funding account at time t, whose value is given by 

B r t ’ = ( 6 ) 

where 

r f ■= r f(v ) = rj l{ y <o} + rfl {y> o } . (7) 

2.2.3 Collateral Account 

Each party posts collateral to mitigate potential losses incurred by the other party in case of default. 
Collateral is in the form of variation margins, i.e. cash transfers are made by the out-of-the money 
party to cover unfavourable moves in the value of the underlying transaction. On {77 A tq > t}, the 
collateral process is given by 

C t := aV(t, St), (8) 

where the collateral level 0 < a < 1 determines the amount of covered exposure. Notice that in our 
setup, if the investor sold the claim to his counterparty, then he is always the collateral taker (Ct < 0), 
while if he purchased the claim from this counterparty he is always the collateral provider (Ct > 0). 
We denote by rj the rate on the collateral amount received by the hedger if he has posted collateral, 
and by r~ the rate paid by the hedger if he has received collateral. The collateral account is defined 
by 

Bj ■= Bj(C) = e /o rc(a)ds , 

where 

r c (x) = r+l {x>0 } + rjt {x< 0} . 

We use to denote the number of shares of collateral account held by the hedger at time t 
holds that 

r t Bj = -C t 


. It then 

(9) 
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2.3 Close-out Value of Transaction 

Since both the hedger and his counterparty can default, the transaction can terminate prematurely. 
In the event of a default, the mark-to-market value of the transaction equals the residual value, 
after mitigating on-default related losses with the available collateral. Let x + := max(x,0), and 
x~ := max(0, — x), be respectively the positive and negative parts of a real number x. Denote by 
0 T {y) the value of the transaction which must be replicated by the hedger at the earlier of the 
two: the hedger’s, or the counterparty’s default time (if positive the hedger owes this amount to the 
counterparty). This is given by 

o T (v) =e T (c,v) (io) 

:= V(r, S T ) + t {Tc<Tl} L c {V(r, S T ) - C T .)~ - 1 {ti<Tc} L 7 (C(t, S t ) - C T .) + 

Here 0 < Lj < 1 and 0 < Lc < 1 are constant loss rates. We refer to Section 3.3 of Bichuch et al. 
(2015) for a detailed explanation. 

2.4 Total Valuation Adjustment 

We review the definition of XVA given in Bichuch et al. (2015). Denote by Vt(<p) the wealth process 
of the hedger associated with the strategy p, where p := (&,£/,£* t > 0). We recall from above 
that C is the number of shares of the funding account. We use £ to denote the number of stock shares. 
Moreover, ^ and £ c denote the number of shares of the bonds underwritten by the trader, and by his 
counterparty, respectively. The wealth process V(p) is given by the following expression 

Vt(<p) := StSt + ZlPt 1 + ttP? + tlBV + - r t B r t c . ( 11 ) 

Bichuch et al. (201 ) show that, under the valuation measure Q, the dynamics of the wealth process 
is given by 


dV t =(r+(gB r t t ) + -rj(Z{B r t t ) + (r D - r~)^ t S t ) + - (r D - r+) (£ t S t ) + r D ^ t P* + r D gP t c ) dt 
~ ( c t) + - r~ (( C t )~) dt + i t aS t dW^ - tf_Pf_ dwl' Q - tf_P t c _dw?' Q 

(12) 

Setting 

Zt = UtaSt, Z\ = -e t -PL Zf = -&_Pf_, 

the dynamics (12) may be rewritten as 

dV t = ( rf(V t + Z{ + Zf - C t ) + - rj (V t + Z\ + Zf - C t )~ 

+ (r D ~ r~)^(Z t ) + - (r D - r+)^(Z t ) - - r D Z{ - r D Z f c - (r+(C t ) + - r~ (Ct)“) 

+ Z t dW + Z{ dwt’ Q + Zf dwf’ Q (13) 

Next, we distinguish between V) + and — V t ~. We use V) + to describe the wealth process of the hedger 
when he replicates the claim &(St) (hence hedging the position after selling the claim with terminal 
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payoff &(St))- On the other hand —V t ~ describes the wealth process when replicating the claim 
— &(St) (hence hedging the position after buying the claim with terminal payoff ^(St))- To this 
purpose, define 


f + (t, v, z, z 1 , z c ; V") = —(r~f (v + z 1 + z c — aVt) + — (y + z 1 + z c — aVt) 

+ (ro ~ r~)-z + - (rn ~ r+)—z~ - rcz 1 - roZ C 
<7 a 

+ r+(aV t ) + - r~{aVt )~) 

f~{t,v,z,^ ,z c -,V) = -f + {t,-v,-z,-^,-z c \-Vt) 


(14) 

(15) 


where the driver depends on the market evaluation process (Vt) (via the collateral (C*)). In particular 
f ± : 14 x [0, T] x M 4 , (oj, t, v, z, z 1 , z c ) eA f ± (t, v, z, z 1 , z c \ Vt(uj)) are drivers of the BSDEs. Moreover, 
define V + , V~ as solutions of the BSDEs 


~dV+ = f + (t, Vt, Z+, Zl’+, Z?' + ; V) dt - Z+ dW? - Z? + dm 1 ® - Z?’ + dm?® 
V+ = 0 T (V)t {T<T} + 4>(S’ t ) 11 {t=t} 


f ~dVf =f-(t,V t -,Zr,Z?-,Z?’--V)dt-ZidW?-Z?-dm I ®~Z?’-dm?^ 
\ v~ = e T (v)t {T<T} + $(S’ T )ii {r=T} 

We have the following main theorem 


(16) 


(17) 


Theorem 2.1 (Bichuch et al. (2015)). Assume that 

r? < r? < r ~, r? < rj , r? V ro < r 1 + hq, rj V < r c + h 1 ?, (18) 

and 

r c v r c < r J < ( r 1 + hq) A ( r c + hJ?). (19) 

If Vq < Vt then there exist prices tt sup and n m t ir in f < ir sup , (called hedger's upper and lower 
arbitrage price) for the claim $(5r) such that all prices in the closed interval [ir tnp , tt sup ] are free of 
hedger’s arbitrage. In particular, we have that ir sup = Vt and ir m f = Vq . 

The total valuation adjustment XVA is defined as the amount that needs to be added to the Black- 
Scholes price (exclusive of funding costs) to get the actual price (inclusive of funding costs). This is 
asymmetric for sell- and buy-prices. Bichuch et al. (2015) define the XVAs as follows. 

Definition 2.2. The seller's XVA is the G -adapted stochastic process (XVAf el1 ) defined by 

XVAf 1 := Vt - V(t, S t ) 


while the buyer’s XVA is defined as 

XVA b t uy := Vf -V(t,S t ). 

XVA sel1 corresponds to the total costs that the hedger incurs when replicating the payoff of a 
claim he sold, whereas XVA buy corresponds to the total costs that he incurs when replicating the 
payoff of a claim he purchased. We note that the difference of the XVAs also describes the width of 
the no-arbitrage interval, as 

XVA s 0 eU - XVA b t y = Vt - Vq . 
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3 PDE representations of XVA 


This section derives the PDE representations corresponding to the master BSDEs (16) and (17). For 
brevity, we will state everything with regards to the upper no-arbitrage price V + , while noting that 
the treatment for V~ is identical. For notational convenience, we will drop the plus superscript, and 
refer to V + , Z + , Z I,+ , Z c,+ and f + simply as V, Z, Z 1 , Z c and /. 

Remark 3.1. In the sequel, we write the solution of the BSDE (16) on {t < r} as v(t, s,w : ,w c ), 
where v(t, St, zu/’®, rof’ 1 ®) = V)ll{ T>t }. The existence of such a measurable function v, i.e. the fact that 
V is Markovian, is shown in Proposition f. 1.1 in Delong (2013). 

To simplify notation, we will not specify the argument of the latter functions. Then we can rewrite 
the BSDE in (16) as 

-dV t = fit, V t , Z t , Z\, Z°- V t )dt - Z t dW® - Z{ dw J t ’ Q - Zfdm?’® (20) 

= f(t, Vt, Z t , Z\, Z t c ; Vt)dt - ZtdW® - f Z{N^idt, dr) - / Z^N c ^{dt, dr), 

J K J M 

V T = 0 T (V) ll {r<T} + T(5 t ) 11 {t=t} (21) 

where j E {/, C}, are the compensated Poisson random measures such that 

wf Q = f [ N^{ds,dr), j E {I,C}, 

Jo Jr 

on [0, r] (we also refer to the proof of Theorem A.2 in Bichuch et al. (2015) for technical details). Let 
v(t, s) be the price of the claim at time t conditioned on St = s, i.e. v(t, St) = Vt. 

Theorem 3.2. Under the no-arbitrage conditions in (18), v is a viscosity solution of the following 
PDE: 

-v t - ^ h < f(0j(v(t,s))-v(t,s,w I ,w c )-v j )-r D sv s -^-<T 2 s 2 v ss (22) 

- fit, V^SVsit^^ 1 ,w c ),6 I {v{t, s)) - v(t, s, w 1 , w c ), 0c(v(t, s)) - vit,s,w J ,w c )-,vit, s)) = 0, 
v(T,s,-,-) = $is). (23) 

Here, we have used the notation v t = jtD, i E {I, C}, and with slight abuse of notation 0i, i E {I, C}, 
corresponds to the price of the contract at default of either party, similar to the way it was defined in 
(10). Specifically, 


Odv) := v + Lciil - a)v) , 

Oiiv) := v — LjU 1 — a)v) + . 

Additionally, v is the unique viscosity solution of the PDE (22) - (23) satisfying the growth condition 
lim |v(-, e x , ■, -)| e _clog2 l a: l = 0, c > 0. 

|x|—>-oo 

Proof. The following proof is an extension of the results from Theorem 4.2.2 in Delong (2013) along 
two directions. The first is that, by contrast with Delong, we use the local property (as opposed to the 
global property) of the viscosity solution (cf. e,g, (Shrevc and Soner, 1994, Section 7)). The second is 
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that we account for the occurrence of defaults, thus we need to extend the argument of Delong (2013) 
to the case when the terminal time is a stopping time, rather than constant. Concretely, we need to 
consider the solution on the set {t < r}. Let £ C 1,2,2,2 ([0, T) xR>o xMxR) be a smooth function, such 
that 4> > v and 4>(to, so, Wq,Wq) = v(to, so, Wq, Wq) for some fixed (to, so,Wq, Wq). Instead of working 
with cj) directly, we will work with a local approximation. Let h > 0 be small enough, such that 
to + h < T. As the default intensities hf, are constant, it follows that Q[r > to] and Q[r > to + h] 
are both strictly positive. Additionally, w.l.o.g. we may assume that cj) is bounded. Specifically, let 
{4>n}n> 1 be a subsequence such that 4> n £ C^’ 2 ’ 2 ’ 2 ([0, T\ x M>o x R x K) is also bounded, satisfying 
4> n A n V (— n) = (j> A n V (— n), and such that it converges to (j> together with its derivatives uniformly 
on compacts. For some no > \v(to, so, Wq, ■u^q 7 ) |, we have that cf> no > v locally. Let U(to, sq, Wq, Wq) be 
an open neighborhood of the point (to, so, Wq, Wq) such that 4> no > v and (by possibly decreasing h ) 
such that (t, so, Wq, Wq) £ U(to, so, w O’ w o) i°r any t G [to, to + h\. Next, choose 

ni>n 0 + sup {|</>(t, s,w J , w c )\, \8c(v(t, s))|, \9j(v(t, s)) |}. 

(t ,S,W J , W C ) GU (to ,so ,Wq ,Wq ) 


This allows us to make the further assumption that (to,so,WQ + 1,Wq) £ U(to, so, Wq, Wq) and 
(to, so, Wo,Wq + 1) £ U(to, so, Wq, Wq) by possibly extending the open neighborhood, because we have 
that (f> ni +i(to,so,WQ + 1 ,wfi) > 0j(v(t, s)), and (f> ni +i(to, s 0 , Wq, + 1) > 9c(v(t,s )). We now have 
that <j) ni +i (t, St, rof’ 1 ®) >Vt,t£ [to, to + ti], even if the default happens during this time interval. 
In the sequel, we will use <j> ni +i in place of cj), but we will abuse notation and continue referring to it 
as cj). 

For to < t < to + h define (V, Z, Z 1 , Z c ) to be a solution to the following BSDE: 


btAril{ r >t} 0 + h) At, S(t 0 -\-h)/\T, W (t 0 +h)Ar’ W (t 0 +h)AT^{T>t} 

r (t 0 +h)Ar _ _ _ „ 

f(r, V r , Z r , zl, Z r c ; V r )dr - 


+ 


/ tAr 


-LQ ^ 

r(t o~\-ti)A- 
It At 


z r dw; - ]T 


r(t 0 +h) At 
It At 


Z^N^(dr, dp). 


This BSDE has a unique solution by Theorem A.2 in lichuch et al. (2015), so V is well defined. By 
the comparison Theorem A.3 in Bichuch et al (2015), we obtain that VtA T ^{ T >t} Z ktArll{r>t}- 
For convenience, we define the operator 

Cu(t, s, w 1 , w c ) = r£)Su s + / ( u(t, s, w 1 + 1 ,w c ) — u(t, s, w 1 , w c ) — ui(t, s, w 1 ,w c ))u I (dz) 

Jr 

a 2 g 2 r 

H-— U ss + / (u(t,s,w : ,w c + 1) — it(f,s,u/,rc c ) — uc(t, s,w : ,w c ))u c (dz), 

2 Jr 

where the Levy measure i/J(dx) = hJSi(dx), j £ {I, C}, with di being the Dirac measure concentrated 
at 1. Additionally, let 

0(t, s, w 1 , w c ) = 4>t(t, s, w 1 , w c ) + Ccj)(t, s, w 1 , w c ), 

T 1 (t, s, w 1 , w c ) = 4>(t, s, w 1 + 1, w c ) — cj)(t, s, w 1 ,w c ), 

T c (t, s, w 1 , w c ) = 4>(t, s, w 1 , w c + 1) — 4>(t, s, w 1 ,w c ), 

V t = V t - cj)(t, S t , ml'^Wt®), 

Z t = Z t - aS t <f>s(t, S t , w 1 ^, Wt3), 

Z 3 t = Zl - P(t, S t , m?l Q ), j £ {I, C}. 
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By Ito’s formula, we have that for to < t < (to + h ) 

i,® „c, 

~tA 

r(t 0 +h)/\T 
It At 


W A + S tAr, ^Ar> W tAr )l{r>t} = <X(*0 + h) A T, 5 (to+ft ) A T, ro (i^ + ft )A r> W ^+h)Ar) %>*} 

c /,Q C,Q W r{to+K)AT 

ky(r,b r ,m r _,w r _)dr — a / 

•Aa 


iAr 


S r (j) s (r, SriVVr’S, W E 


- E 

ie{/,c} 


r(t 0 +h)AT 
It At 


[ f J (r 

Jl 


r, SV, ,w^A)N :i '' l ^(dr, dp). 


It follows that 
BtArl{r>t} = 


r(t 0 +h)AT . _ 

I (@(r,Sr, 

tAT V 

f(t 0 +h)Ar = 

Z r dW.f - J2 


wl’S, w°®) + f(r, Vr + <t>,2 r + aS r (f> s , E + F, Z? + F 7 ; V r ))dr 

•(to+h)AT 


' ZAr 


ie{/,C} 


It At 


Z ] r N M (dr, dp). 


(24) 


Now, assume by contradiction that 4> violates the subsolution property, i.e. there exists e > 0 such 
that 


- cj)t(to, s 0 , Wq, Wq) - £4>(t 0 , s 0 , Wq, Wq) 

- f(to, crs 0 (j)s(to, s 0 , WQ, wu),!: 1 (t, 0 , s 0 , w^, Wo),f c (t 0 , s 0 ,Wo,Wo);v(to, s 0 )) > e. 

By continuity, and by possibly reducing h and the open neighborhood of (to, so, Wq, Wq ), we may 
assume that it is also true inside the neighborhood Id (to, sq, Wq, Wq) of (to, so,Wq,Wq), while still 
assuming that (t, so, Wq, Wq) £ Id (to, so, Wq,Wq) for any t £ [to, to + h]. 

Next, define 


h : = 


r(to+h)A7 


Jt 0 AT 


\H( t , St, wl'^wf'^dt 


where 


^>(t, s, w 1 , w c ) :=(j>t(t, s, w 1 , w c ) + Ccj)(t, s, w 1 , w c ) 

+ f(t, 4>, crs<j) s (t, s, w 1 , w c ), f \t, s, w 1 , w c ), f c (t, s, w 1 , w c )\v(t, s)). 

Since / is Lipschitz, we can show in a similar way as in dong (2013) that for to < t < to + h it holds 
that 


\H(t, s, w 1 , w c ) <(7(1+ s 2 ), 


V t 


tAT 


1 


{r>t} 


< Chi 

(t 0 +h)AT - 2 

z. 


dr + ^ 


r(t 0 +h)AT r 

_ 

2 . 

JtoAr J R 

zi 

zd (dp) dr 


< Chi, 


/i°Ar ^ ^ 

< (to + h) At] < Ch, 
where we have introduced another stopping time 

ti := inf {t > t 0 : (t, S t , zcrf’®, rof ,<Q ) 0 U(t 0 , s 0 , Wq, Wq)}- 


(25) 

(26) 
(27) 
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From q\ := Q[r > to] > 0, inequality (27) and the fact that the default intensities are 

constant, so that Q[t £ (to, to + h)\ < Ch, it also follows that 

Q[{ri > (to + h) A r} n {t > to + h}] 

> 1 - Q[t < t 0 ] - Q[n < (t 0 + h) A r] - Q[t £ (t 0 , t 0 + /i)] > qi - Ch, 
and similarly that 


Q[{n < (to + h) A U {r £ (t 0 , to + h)\\ < Ch. 

Here and throughout this proof, we will abuse notation and use C > 0 to denote a generic constant, 
which may be different in each inequality below. It follows that 


Ih = 


h 


+ 


r(to+h)A7 

Jto At 

r{t 0 +h)An 

Jto At 


.y r /. Q 7 x ^7 J-Tl 

^(T St, )utH{ Tl >(t o +/i)AT}n{T0(t o ,to+h)} 

.y./j. q _CiQAjj.ii 

^(T St j i^t / ( 7^ii{ri<(to+/i)Ar}U{re(io,<o+7*)} 


< — e(Q>[{n 

> (t 0 + h) A r} n {t > to + /i}] 




■ rto+h 

/ ldtl{ r i<(i 0 -)_/ l )Ar}U{re(to,to+7i)} 
■ J to 


' r(t 0 +h)AT 

/ ^ 2 (i,St,OT*’ Q ,tuf’ Q )cft 

.JtoAT 

< —eQ [{ri > (to + h) A r} n {r > t 0 +}] 




+ 


■ Mto+h) At 
HqAt 


dt 2 (t, St, vjf'^)dt 


< -e(?i - Ch) + 


\ 


1 + EQ 

sup Sf 


Ue[toAr,(to+h)Ar] J 


This shows that for sufficiently small h > 0, we have that Ih < — 

From the comparison Theorem A.3 in Jichuch el al. (201. ) we get that Vt 0 ArH{r>t 0 } — 0. Using 
this along with (24), we conclude that 


e 

- < 
2 “ 


^^toArl{T>t 0 } Ih 


' r(to+h)AT — =_=_ __ 

/ /(r, Uy + 0, z r + uSV^, Zl + f ■ J , Z r c + f c ; U r ) - f(r, 0, aSy</> s , f 7 , f c ; Uy)dr 

Jto At 


< C sup 

r£[to,to+h] 


Vr 


tAt 




{r>r} 


+ 


c 


MV 

< C(h^ + h^), 


r(to+h)A7 
I to At 


Z r 


dr 


+ £ 


^■(io+ft)Ai 
.JtoAT 


• 2 


Zy i/J(dp)dr 


where the second inequality follows from the fact that the driver / of our BSDE is Lipschitz along 
with the Holder’s inequality (setting the exponent to 2), while the last inequality follows from the 
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estimates given by the inequalities (25) and (26). This yields a contradiction, as we let h \ 0. Hence 
<p does not violate the subsolution property. 

In order to incorporate the jump condition (21) at the jump time r, note that we can rewrite the 
subsolution property for cf as 

- Mto, so, u’o, Wo) + h j4>j{t Q , so, w'o, w o ) “ r DS(j)s ~ 2 s 2 (j) ss (28) 

- <j>, CrS 0 4>s(to, S 0 , Wq, wfi ), f' \t 0 , S 0 , Wq, wfi),T C (to, S 0 , Wq, wfi); v(t 0 , s 0 )) < 0, 

where we set f\(t, v, z, z 1 , z c ; v(t, s)) := f(t, s, v, z, z 1 , z c , v(t, s)) + lifz 1 + h^z c . Under the no arbi¬ 
trage conditions given by (18), it follows that fi is an increasing function in both arguments z 1 and 
z c . By continuity of v (see e.g. Lemma 4.1.1 of Delong (2013)), it follows that <j>(ri, S TI , + 
l,wf; Q )ll {rj < TcAT} > V Ti 1 {ti<tcAT} = 0j(V Tl )H{r/<rcAT}, and similarly ^>(r c , S TC , w™ + 

1 )i{r c <r I AT} > l{r c <r / AT} = 0c(Vr c ) 1 { Tc<TjA t} • Together with (28) it follows that 

- <k(to,s Q ,Wo,w$) + h’f&jito, s 0 , Wq,Wq) — r D s(f> s — -s 2 cr 2 (f) ss 

je{i,C} 

- hf{0j(v) - (f>{t 0 , s 0 , Wq , Wq" )) - h2(0 C (v) - 0(t o ,s o , Wo, Wo)) 

- f(to, </>, o-s 0 (/) s (to, So, Wq,Wq), 0/(v) - ^{to,s o ,Wo,Wo),0c{v) - <j>(to,So,Wo,Wo)\v(to,So)) < o. 

This shows that 0 is a subsolution of the PDE (22) as claimed. 

Finally, the uniqueness result, follows from the uniqueness of the solution to the BSDE, which in 
turn follows from the comparison Theorem A.3 in (2015). □ 

Remark 3.3. Since we are only concerned with Vj before any default occurs, there is no need to keep 
track of the martingale terms ’s. These are only needed to realize that a default has happened. 
Otherwise, zu 3 /^ = —hjt. It then follows that v reduces to a function of only two variables, so that we 
can simply define v(t,St) = I n other words, v(t,s) = v(t, s, —loft, —h^t). In this case, the 

PDE (22)-(23) becomes 

-v t + (/*/ + h$)v(t, s) - r D sv s - ^ c t 2 s 2 v ss (29) 

- f(t, v, asv s (t, s),9j(v(t, s )) - v(t, s)d c (v(t, s)) - v(t, s)\v(t, s)) = ^ hj0j(v{t, s)), 

je{i,C} 

v(T, s) = <&(s). (30) 


Remark 3.4. Additionally, we will also employ the standard change of variables x = logs, so that 
w(t,x) = v(t,e x ) and w(t,x) = v(t,e x ). Note that this change of variables allows us to get rid of 
boundary conditions at S = 0. Then, the PDE (29) together with the boundary condition (30) becomes 


-w t - (r D - y) w x - \wr 2 w xx + (hf + h^)w 
- f(t,w,aw x (t,x),9i(w(t,x)) - w(t,x),0 c {w(t,x)) - w(t,x)\w{t,x)) 

w{T,x) = ${e x ). 


(31) 

ie{i,c} 

(32) 
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Hence, we can express the whole pricing problem as a Cauchy problem for a two-dimensional system 
of semilinear PDEs: 


—wt + C l w = f(t, w, c 7w x , Oj(w) — w, 0c(w) — w: w) 
—wt + C 2 w = 0 
w(T,x ) = w(T,x ) = < &(e a: ), 


where the differential operators are defined by 


■= -Ud ~ -tt)9 x ~ %~d xx + (hf + h%) ■ - 

£ 2 := -(td ~ y)&r - ydxz 


It turns out that the PDE (31)-(32) (resp. (29)-(30)) not only has a unique viscosity solution, 
but this solution is also a classical solution if we assume that $ and < h / (where defined) have at most 
polynomial growth, i.e., |$(s)| < C( 1 + s n ), |^'(s)| < C( 1 + s n ), s G M>o The classical argument of 
Theorem 20.2.1 in 'annon (19 ) assumes a bounded terminal condition, but we can employ a change 

of variables and divide by (1 + s 2n ) to utilize their framework: 


= (l + e 2nx )w{t,x), ti(t,x) = (l + e 2nx )w{t,x), $(x) = (1 + e 2nx )<f>(e x ). (33) 

In this case, our PDE (29)-(30) becomes: 


1 2 = 

-w t - -c r w xx + 2 na 


s (2n—l)x 




‘ X + e 2nx w * ~ rDW *) 


(34) 


+ [hY + h^ + 2 nr D 


D (2n—l)x 


D (2n—2)x 


2 nx 


l + e 

,(2n— l)x 


+ na“ 


l + e 


2 nx 


i l - 2n - ^ 


w 


( ^{Zn—l)x _ _ ^ \ 

t, w, aw x - 2 na i 2nx w,e!(w(t,x)) - w, 9 c {w(t, x)) - w;w(t,x)j = hj9j{w{t,x)), 




w(T, x) = 3>(x). 


(35) 


We can then prove the following 

Proposition 3.5. Assume that is piecewise continuously differentiable and $ as well as <h / (where 
defined) have at most polynomial growth, i.e., |<J>(s)| < C( 1 + s n ), |<P'(s)| < (7(1 + s n ), s G M>o- Then 
the PDE (31) with terminal condition (32) has a classical solution. 

Proof. Assume first that $ is continuosly differentiable. Using the transformation (33) given above, it 
is then suffices to prove that (34)-(35) has a classical solution. The above transformation guarantees 
that both and $>' are bounded. Then the existence of a smooth (and bounded) solution to (34)-(35) 
follows from Theorem 20.2.1 in 'annon (195 ). In the case that $ is only piecewise smooth, the 
original proof can be modified following a similar procedure to Jouini and Kallal ( I! ). Hence, using 

the change of variables (33), we conclude that there exists a classical solution to the PDE (31) with 
terminal condition (32). □ 
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Combining Proposition 3.5 with the uniqueness result from Theorem 3.2 allows us to conclude that 
there exists a unique classical solution to the PDE with polynomial growth. Moreover, the hedging 
strategies can be easily found by an application of Theorem A.l. Using the relations in (13) we find 
that on the set {t < r} 


& = = vs(t,S t ), 

abt 

j _Z[ = v(t, S t ) - OjtvfaSt)) 

p{ e -(r D +hJ)(T-t ) ’ J 

where we have used that the bond price P( = e^ rD+h i ^ on tj > t. 
between the measures P and Q given in (3). 


(36) 

by virtue of (2) and the relations 


4 Numerical Analysis 

We perform a comparative statics analysis to analyze the dependence of XVA and portfolio replicating 
strategies on funding rates, default intensities, and collateral levels. We consider the relative XVA, 
i.e. express the adjustment as a percentage of the price V) of the claim, given by 1 '*A* , where 

Vt = Vp depending on whether we are considering buyer’s or seller’s XVA. The claim is chosen to be a 
European-style call option on the stock security, i.e. 4>(x) = (x — K) + ■ We consider one at-the-money 
option, with So = K = 1 maturing at T = 1. In order to focus on the impact of funding costs (which 
in practice is the most relevant) and separate it from additional contributions to the XVA coming 
from asymmetries in collateral and repo rates, we set = r~ = 0.05, and = r~ = 0.01. As the 
derivatives contract does not only specify the price of the option but also the levels of collateralization 
of the deal, the no-arbitrage region appears as a (two-dimensional) band in XVA and a rather than 
as a (one-dinrensional) interval in XVA only. 

We also use the following other benchmark parameters: <r = 0.2, = 0.05, rj = 0.08, rp = 0.01, 

r 1 = 0.03, r c = 0.04, hf = 0.2, h® = 0.15, Lj = Lq = 0.5, and a = 0.9. We compute the numerical 
solution of the PDE using a finite difference Crank-Nicholson scheme. 

The main finding of our analysis are discussed in the sequel: 

Higher funding rates increase the width of the no-arbitrage band. Figure 1 displays the no¬ 
arbitrage band whose width is increasing in the funding rate rj. As a gets higher, the band noticeably 
shrinks reaching its minimum around a = 80% before widening again. Notice that buyer’s and seller’s 
XVA do not have a symmetric behavior. This can be better understood by analyzing the dependence 
of the band on the collateral level a in Figure 1. If a is not too high (a < 0.5), the widening of the 
no-arbitrage band with respect to the funding rate rj is due to decreasing buyer’s XVA. On the other 
hand, if a is high the buyer’s XVA is insensitive to changes in rj whereas the seller’s XVA increases 
with rj, contributing to widen the no-arbitrage band. This is further supported by the numerical 
values reported in Table 1. When a < 0.5, the position in the funding account for the seller’s XVA is 
long and the same regardless of the funding rate rj. On the other hand, the size of the long position 
for the buyer’s XVA increases in rj. In presence of full collateralization, i.e. a = 1, the situation 
reverses. The position in the funding account for the buyer’s XVA is short and stays constant with 
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respect to rj . Vice versa, for the seller’s XVA the size of the short position increases in size with 
respect to rj. 

If a is high, the trader will have to post more collateral and consequently reduce the cash resources 
for his replicating strategy. He will then have to borrow more from the funding desk, resulting in higher 
funding costs. This drives up both the seller’s XVA and the number of shares of stocks and bonds 
needed for the replication strategy. 



a a 



Figure 1: Top left: Buyer’s and seller’s XVA as a function of a for different rj . The seller’s lies above 
the buyer’s XVA and the same line style is used for both. Top right: Number of stock shares in the 
replication strategy. Bottom left: Number of trader bond shares in the replication strategy. Bottom 
right: Number of counterparty bond shares in the replication strategy. We plot the strategies for the 
portfolio replicating the seller’s XVA. 


Higher collateralization increases portfolio holdings. As the collateral level a increases, the 
seller’s XVA increases. This happens because the value of the closeout position becomes higher as 
it can be directly seen from Eq. (10) (notice that V > 0 because we are considering a short call 
option position). The trader would then need to construct a portfolio replicating a larger position, 
hence he must take more risk. He achieves this by increasing the number of shares of stock and 
bond underwritten by the counterparty. Moreover, higher collateralization levels reduce the size of 
the downward negative jump to the closeout value occurring when the trader defaults. Consequently, 
the trader needs to purchase a smaller amount of his bonds to replicate this position as a increases. 
This behavior is confirmed from the plot in Figure 2. 
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a 

r f 

Seller’s XVA: funding account ($) 

Buyer’s XVA: funding account ($) 

0 

0.08 

0.0039 

0.0403 

0 

0.2 

0.0039 

0.0447 

0.25 

0.08 

0.0249 

0.0257 

0.25 

0.2 

0.0249 

0.0287 

0.75 

0.08 

-0.0037 

-0.0036 

0.75 

0.2 

-0.0038 

-0.0032 

1 

0.08 

-0.0182 

-0.018 

1 

0.2 

-0.0193 

-0.018 


Table 1: The columns give the dollar position in the funding account corresponding to the replicating 
strategies of seller’s XVA and buyer’s XVA. 



Figure 2: Top left: Buyer’s and seller’s XVA as a function of a for different h®. Top right: Number of 
stock shares in the replication strategy. Bottom left: Number of trader bond shares in the replication 
strategy. Bottom right: Number of counterparty bond shares in the replication strategy. 


The width of the no-arbitrage band is insensitive to counterparty’s default intensity. 

Figure 2 shows that both seller’s and buyers’s XVA decrease, if the counterparty’s default intensity 
h® increases. When a is low, the two quantities drop by nearly the same amount and the width of 
the no-arbitrage band is unaffected. As a gets larger, the seller’s XVA decreases faster relative to the 
buyer’s XVA and the two quantities almost coincide when a = 1. 
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r f 

Seller’s XVA: funding ($) 

Buyer’s XVA: funding account ($) 

0.08 

-0.0124 

-0.0123 

0.1 

-0.0125 

-0.0122 

0.15 

-0.0127 

-0.0122 

0.2 

-0.013 

-0.0122 


Table 2: The columns give the dollar position in the funding account corresponding to the replicating 
strategies of seller’s XVA and buyer’s XVA. We set h® = 0.15. 


Consistently with Figure 2, Figure 3 shows that the seller’s XVA decreases when the default 
intensity of the counterparty h® increases. This can be understood as follows. Using the relations 
between the default intensities under the probability measures P and Q given in Eq. (3), we can see 
that increasing h® is equivalent to increasing the bond rate r c while keeping the default probability 
under P constant. Hence, the trader would earn higher premium from his long position in counterparty 
bonds (see also bottom panels of Figure 3). Such a gain dominates over the funding costs incurred 
when replicating a larger closeout position (Eq. (10) indicates that the closeout payment increases to 
the risk-free payoff V as /t<S increases). Altogether, this means that the funding costs of the investor 
would be reduced as /i® increases. 

Hence, the dependence of XVA on counterparty’s default intensity contrasts with its sensitivity 
to collateral levels numerically illustrated in Figure 1. This is because higher default risk of the 
counterparty also means higher return on the bond underwritten by the counterparty, whereas higher 
a only means that a larger value of the closeout payment needs to be replicated. 

5 Conclusions 

We have developed a rigorous analysis of the semilinear PDE associated with the BSDE characterizing 
the price process of a portfolio replicating a European option, when funding, collateral and closeout 
costs are taken into account. We have shown the existence and uniqueness of a classical solution to 
the PDE under mild assumptions on the coefficients. Using this result, we have conducted a thorough 
numerical study analyzing the sensitivity of XVA and of the claim’s replication strategy to collateral 
levels, default risk and rates asymmetries. Our findings support the introduction of centralized XVA 
desks to manage and hedge all costs related to over-the-counter transactions. It shows that funding 
costs originating from the different trading components cannot be easily separated and hence attributed 
to different business units (CVA, DVA and FVA desks) because they are highly interdependent. 

A PDEs and Replication Strategies 

Recall the measurable function v from Remark 3.3 that was defined by v(t,St) = Vtll{ T >n. The 
following theorem shows how the function v can be used to compute the hedging strategies. Noting 
that the law of St is absolutely continuous, the proof of the theorem becomes analogous to the proof of 
Theorem 4.1.4 in iclong (2013). Hence, we omit it here and only give the statement of the theorem. 
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Figure 3: Top left: Seller’s XVA as a function of h® for different a. Top right: Number of stock 
shares in the replication strategy. Bottom left: Number of trader bond shares in the replication strat¬ 
egy. Bottom right: Number of counterparty bond shares in the replication strategy. The replicating 
portfolio refers to the seller’s XVA. 


Theorem A.l. Consider the data (/, 9 t (V(t, S t )),V(T, St)) for the BSDE given by (16). Addition¬ 
ally, let the function v(t,St ) = V)ll{ T> t} be defined as in Remark 3.3. Then, on the set {f < t} we 
have that 

Z t = aS t vs(t, S t ), 

zi = e j (v(t,s t ))-v(t,s t ). 
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